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We have performed high resolution angle resolved photoemission (ARPES) studies on electron 
doped cuprate superconductors Sm2-xCe x Cu04 (x=0.10, 0.15, 0.18), Nd2-ccCe x Cu04 (x=0.15) 
and Eu2— ccCe^CuC^ (x=0.15). Imaginary parts of the electron removal self energy show step- like 
features due to an electron-bosonic mode coupling. The step-like feature is seen along both nodal 
and anti-nodal directions but at energies of 50 and 70 meV, respectively, independent of the doping 
and rare earth element. Such energy scales can be understood as being due to preferential coupling 
to half- and full-breathing mode phonons, revealing the phononic origin of the kink structures. 
Estimated electron-phonon coupling constant A from the self energy is roughly independent of the 
doping and momentum. The isotropic nature of A is discussed in comparison with the hole doped 
case where a strong anisotropy exists. 

PACS numbers: 74.25.Jb, 74.72.-h, 79.60.-i 



Electron-bosonic mode coupling in solids manifests it- 
self as a slope change in the dispersion or a "kink" in 
angle resolved photoemission (ARPES) dispersions. 1 
ARPES studies have been extensively performed to 
obtain information on the electron-bosonic mode cou- 
pling (EBC) 2 '-^^* 7 - in high temperature superconduc- 
tors (HTSCs) to find out what mediates the electron pair- 
ing. In spite of the extensive studies, controversy still ex- 
ists on what causes the strong kink structure in ARPES 
spectra. The controversy can be rooted in the fact that 
there exists two bosonic modes with similar energy scales 
in hole doped HTSCs that may couple to quasi-particles, 
i.e., phonons and magnetic resonance modes. Optical 
phonons have energies between 40 and 90 meV 8,9 while 
the magnetic resonance mode from neutron experiments 
shows an energy scale of 40 me V 10 i n . Various devised 
ARPES experiments have been performed to pin point 
which mode causes the kink effects. However, the con- 
troversy will continue because the two energy scales are 
similar. 

The magnetic mode in electron doped HTSCs mean- 
while has been discovered only recentl y 12 ! 13 . While the 
energy scales of the optical phonons in the Cu02 planes 
are similar between electron and hole doped cuprates, 
the magnetic mode in electron doped cuprates is found 
to be much smaller (~ 10 meV or smaller). Therefore, 
it should be easier to discern the two effects in electron 
doped HTSCs. In addition, recent progress in the crystal 
growth technique provides high quality samples suitable 
for EBC studies 14 . This recent progress gives us an im- 
portant opportunity to study the kink effects in electron 
doped HTSCs. In fact, there is a recent report on this 
issue from the scanning tunnelling spectroscopy 15 . Moti- 
vated by this issue, we performed ARPES experiments 
on various electron doped systems of Sn^-^Ce^CuCU 
(SCCO), Nds-^Ce^CuCU (NCCO) and Vu 2 - x Ce x CuO A 
(ECCO) with specific aim on the EBC studies. Our re- 



sults show a clear evidence for electron-phonon coupling 
(EPC). Surprisingly, the coupling is isotropic in a strong 
contrast to the hole doped case. This fact should help 
to advance our comprehensive understanding of HTSCs 
over the entire phase diagram. 

SCCO 0=0.10, 0.15 and 0.18), NCCO (x=0.15) 
and ECCO (x=0.15) single crystals were grown by the 
travelling- solvent floating-zone method. Relative Ce con- 
centration was determined by Ce core-level photoemis- 
sion and found to be consistent with the nominal val- 
ues. NCCO, SCCO and ECCO were reduced by anneal- 
ing in N2 for 10 hours at 960, 900, 850 C, respectively, 
and then in oxygen for 20 hours at 500 C to induce su- 
perconductivity. T c was determined to be 0, 17, 9, 23 
and K by magnetic susceptibility measurements for 
SCCO(x=0.10), SCCO(x=0.15), SCCO(x=0.18), NCCO 
and ECCO samples, respectively. ARPES experiments 
were performed at beamline 5-4 of the Stanford Syn- 
chrotron Radiation Laboratory using 16.5 eV photons 
with an energy resolution of 14 meV. Samples were 
cleaved in situ and laser aligned. The chamber pressure 
was better than 4 x 10 -11 torr and the temperature was 
kept at 15 K. 

Fig. 1 (a)-(d) shows ARPES data along the nodal di- 
rection from the SCCO samples with £=0.15 and 0.18. 
x=0.1 data is not shown because the feature is too broad 
to discuss the EBC. The raw data and the momentum 
distribution curve (MDC) dispersions show kink-like fea- 
tures at around 50 meV (marked by the arrows), strongly 
suggesting an EBC at that energy. However, it was shown 
that electronic structures of electron doped HTSCs show 
band folding effects due to possible \[2 x y/2 orderin g 17 ! 20 . 
In fact, x=0.15 data has a nodal gap, as shown in the in- 
set of Fig. 1 (b), stemming from a band folding effect 14 . 
In such case, it is difficult to judge if the bend in the dis- 
persion is a kink or from band folding effects, especially 
for the nodal region. 
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FIG. 1: (Color online) Intensity maps of ARPES data along 
the nodal direction from (a) x=0.15 and (c) 0.18 samples of 
SCCO. The lines in panels are the MDC dispersions while 
the arrows mark the energy where the kink structure appears 
in the dispersion. Panels (b) and (d) plot the EDC stacks 
of the data near the top of the band. The vertical dash- dot 
line marks the kink energy. The inset in panel (b) shows the 
lowest binding energy feature which remains below the Fermi 
level. The inset in panel (d) depicts the direction of the cuts 
in the Brillouin zone. Simulation results on EBC for a folded 
bare band is shown in panel (e) with EDCs in panel (f) in a 
similar fashion to the experimental data. 



To see if one could still observe a kink structure in spite 
of a band folding effect, we plot in Fig. 1 (e) a simulated 
ARPES spectral function A(k,cj) of a folded band with 
a 50 meV Einstein phonon. The bare band was taken 
from the published nodal dispersion of SCCO^ and a 
moderate EPC constant A = 0.5 was used. One finds 
that the kink structure is still visible for a folded band 
with a moderate A. 

We will later provide discussions on the EBC based on 
the self energy analysis. However, one can already find 
some EBC aspects from the energy distribution curve 
(EDC) plots. Shown in Fig. 1 (f) are EDCs from the 
simulation data near the top of the band. Not only can 
one still see a kink-like dispersion, but also observe that 
the peak becomes considerably broad beyond the mode 
energy as scattering is allowed through EBC. Looking at 
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FIG. 2: (Color online) Intensity maps of ARPES data along 
the anti-nodal direction from (a) x=0.15, (b) 0.10 and (c) 0.18 
samples of SCCO. As in Fig. 1, the lines are MDC dispersions 
and the arrows mark the energy where kinks in the dispersion 
exist, (d) EDC stacks of the x=0.15 data in (a). The vertical 
dash- dot line marks the kink energy. 



the experimental EDCs in Fig. 1 (b) and (d), one can 
also see that the peaks quickly broaden beyond the ener- 
gies marked by the vertical dash-dot line. Such behavior 
already suggests the existence of EBC in these materials. 

The kink feature is not only observed along the nodal 
direction but also in the anti-nodal direction. Fig. 2 
shows ARPES data along the (7r,0) to (7r,7r) direction. 
The data show clear Fermi edges unlike the spectra from 
the nodal direction. This is because the near Ej fea- 
tures are from momentum points that are far enough 
from the anti-ferromagnetic Brillouin zone (BZ) bound- 
ary and thus the folding effect is minimal. The MDC 
dispersions again show kinks at the arrow-marked en- 
ergy positions. With minimal folding effects in the re- 
gion, one can attribute the kinks to the EBC. EDCs of 
£=0.15 Fig. 2 (d) also show that sharp peaks at the low 
binding energy quickly become broad beyond a certain 
energy, which is a clear sign of EBC. It is noted that the 
characteristic energy where the kink exists is at around 
70 meV (compared to 50 meV for the nodal direction) 
and is more or less doping independent. 

A natural and important question is whether such 
kink-like feature is an intrinsic property of electron doped 
HTSCs or is specific to SCCO. In that regards, investi- 
gation of other electron doped compounds is important. 
Shown in Fig. 3 are the data from optimally doped 
(£=0.15) NCCO and ECCO. The NCCO nodal cut in 
Fig. 3a shows a dispersion that crosses the Fermi en- 
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FIG. 3: (Color online) Intensity maps of ARPES data, (a) 
Nodal cut for NCCO 0=0.15) and (b) EDCs. Anti-nodal 
cuts for (c) NCCO and (d) ECCO (x=0.15). 



ergy, unlike the SCCO with £=0.15 case. The EDCs 
of the data presented in panel (b) again show sudden 
broadening beyond a characteristic energy. On the other 
hand, ECCO nodal dispersion has a gap of ~ 4QmeV 
and spectral features are broad. Therefore, for the same 
reason as the x=0.10 SCCO case, ECCO nodal data is 
not presented. This trend of the nodal gap increase when 
Nd is replaced by Sm and Eu is consistent with earlier 
observation 16 . For anti- nodal cuts, depicted in Fig. 3 (c) 
and (d) , the dispersive features cross the Fermi level and 
show clear kink-like features near 70 meV as was the case 
for SCCO. The above results tell us that EBC is present 
in all the electron doped HTSCs, independent of the rare 
earth element. 

For a quantitative analysis on EBC, it is useful to look 
at the self energies of the spectral functions. Convention- 
ally, the real part of self energy is obtained by subtracting 
the bare band dispersion from the MDC dispersion and 
imaginary by multiplying the MDC width with the Fermi 
velocity. This method, however, is not applicable to the 
data from these materials, especially for the nodal data, 
because the bare band dispersion is not known due to 
a folding effec t 14 ! 20 . One could think about fitting the 
ECDs, but it is too unreliable due to the highly asym- 
metric lineshape. 

In such cases, one may use a newly developed method 
for determining the electron removal self energy 21 . 
ARPES intensity l(k,u) is proportional to the imaginary 
part of electron removal Green's function or the spec- 
tral function A(k,u). One could recover the full Green's 
function if one knew A(k,w) over the entire energy range 



because the real and imaginary parts of the Green's func- 
tion are related through the Hilbert transform. However, 
this is not possible because ARPES measures the spec- 
tral function only below the Fermi energy (uj=0). In 
that case, one can define electron removal Green's func- 
tion which bares information on the electron dynamics 
only below lj=0. The rest comes naturally and one can 
obtain the imaginary part of the electron removal self en- 
ergy ImS^ without any a priori assumptions. ImS^ is 
then a simple algebraic expression of the ARPES spec- 
tral function and its Hilbert transform 2 ^. The advantage 
of the method is that it can be applied to any line shape 
and represents the "peak width" of a spectral function. 
For the purpose of the discussions given here, ImS^ can 
be regarded as ImS^i. 

Fig. 4a shows ImS of the nodal cuts in Figs. 1 and 3. 
The data in panel (a) show a considerable slope change 
around 50 meV (vertical dashed line). For the data along 
the anti- nodal direction in panel (b), the slope change oc- 
curs around an increased energy of 70 meV. While the 70 
meV kink observed near the anti-nodal direction has been 
observed in earlier studies on NCCO 22 , the kink at 50 
meV along the nodal direction in electron doped HTSC 
is seen for the first time 2 ^. The slope change around 20 
meV is an artifact due to the finite energy resolution. 

Our data show that there is EBC for both nodal and 
anti-nodal directions, at ~50 and ~70 meV, respectively. 
These mode energies are roughly independent of doping 
and rare earth atoms. Then the question is which bosonic 
mode causes the observed kink structures. Recently, it 
was found through neutron experiments that the energy 
of the spin resonance mode in electron doped HTSCs 
is at most 10 me V 12 i 13 . Therefore, the spin resonance 
mode can not explain the kinks at 50 and 70 meV in the 
ARPES spectra from electron doped HTSCs. 

Instead, it can be naturally attributed to the bond 
stretching mode of the copper oxygen plane. Fig. 4(c) 
shows the published phonon dispersions of bond stretch- 
ing modes in electron doped HTSCs 2 ^. We note that the 
energy range of the phonon dispersion is similar to that 
of the kink energy. It is known that electron coupling 
to bond stretching phonons is stronger near the BZ zone 
boundaries. Phase space argument can be used to show 
that an anti-nodal photo-hole couples to full breathing 
phonons (#=[0.5,0.5,0], 60 meV) 25 and experimental re- 
sults show that a nodal photo-hole mostly couples to half 
breathing mode (#=[0.5,0,0], 50 meV)^^. This is very 
much compatible with the observed energy scales. 

As it is clear from the above discussion that the kinks 
in our data are from the EPC, we try to estimate the 
EPC constant. Conventionally, the EPC constant A is 
obtained from the mass re-normalization factor, that is, 
the ratio between the Fermi velocities of the bare and 
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this method does not work when the band dispersion is 
non-linear, which is the case for SCCO. On the other 
hand, one can estimate A from the ImE in a more gen- 
eral way. The step height is already a rough measure 
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FIG. 4: (Color online) (a) Imaginary parts of self energies 
ImE of the nodal data in Fig. 1. Curves are offset by 50 meV 
for clarity. The bottom curves show phonon contribution to 
ImE after a linear background has been subtracted (solid) 
and theoretical ImE with an Einstein phonon (dash-dot), (b) 
ImE for the ant i- nodal direction with 30 meV offset, (c) Ex- 
perimental phonon dispersions of the bond stretching modes 
in NCCO. (d) Estimated EPC constants from ImE in panels 
(a) and (b). 



of EPC strength. To obtain the step height, one has to 
subtract the electron-electron interaction contribution to 
ImE. For the energy window in Fig. 4a, ImE has an 
approximate linear background above the kink energy. 
The solid blue curve is the resulting ImE after the lin- 
ear background has been subtracted, the resulting ImE 
is compared with an ideal ImE from a single Einstein 
phonon at the kink energy (dash-dot line). Even though 
the two are different as there are multiple phonon modes, 
we use the single Einstein model for a rough estimate of 
EPC constant. In that case, the EBC constant is ob- 
tained by a simple equation given in Fig 4a. The results 
are plotted in Fig. 4d. Overall, the value is about A=0.8 



for both nodal and anti-nodal directions and does not 
vary much due to doping and rare earth atom substitu- 
tion. Even though the error bars are as large as 40%, this 
clearly casts a sharp contrast with the hole doped cases 
where EBC constant in the anti-nodal region is as large 
as ten times that of the nodal region^. We can therefore 
say that EPC in electron doped cuprates is (relatively) 
isotropic. 

The fact that EPC universally exists in electron doped 
cuprates and is similar independent of doping and rare 
earth elements while the T c varies from 25 to K suggests 
that the observed EPC in the present work may not play 
the dominant role in determining the superconductivity. 
This fact in turn implies that there must be other impor- 
tant factors in the superconductivity. The clue to this 
question could be found in comparing the present results 
with that of hole doped materials for which T c 's are much 
higher. First of all, contrary to the electron doped ma- 
terials, the EBC A is strongly anisotropic for hole doped 
cuprates. Such anisotropy comes from from the fact that 
a bosonic mode at 40 meV, either B^ phonon 5,25 or 
magnetic resonance mode 3 , preferentially couples to the 
electrons in the anti-nodal region. Electron coupling to 
this bosonic mode is much stronger than the coupling to 
the bond stretching phonon mode 4,25 , and gives strong 
anisotropic EBC in hole doped HTSCs 4 . Secondly, the 
strong EBC in the anti- nodal region of hole doped HTSCs 
is observed only bellow T c 3 while the kink structure still 
exists even in the normal state of electron doped HTSCs 
as can be seen in ECCO. Temperature dependence stud- 
ies on NCCO also show no change in the kink structure 
across the T^ 

These observations point to the fact that there is a 
missing bosonic mode in electron doped HTSCs that cou- 
ples to the electrons in the anti-nodal region and causes 
strong anisotropy in hole doped materials. This dispar- 
ity between the two systems may be related to the fact 
that the energy position of the van Hove singularity is at 
about 350 meV for electron doped while it is very close to 
the Fermi energy for hole doped materials. Independent 
of the origin, the missing mode could be closely related 
to the disparity in T c s between electron- and hole-doped 
HTSCs. 
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